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13  ABSTRACT 


A  standing  wave  acoustic  parametric  source  (SWAPS)  is  proposed  as  a  small,  low -frequency  underwater 
sound  source.  SWAPS  is  a  liquid-filled  cylindrical  tube  which  is  driven  ?.i  one  end  by  a  piston  transducer  and 
terminated  on  the  other  end  b\  a  pressure-release  reflector  to  form  a  resonant  cavity.  The  piston  is  driven  simul¬ 
taneously  at  two  high  frequencies  ^  and  cj2  which  are  at  or  near  resonance  for  the  plane-wave  mode  in  the 
cavity.  The  resulting  large-amplitude  primary  sound  waves  mi.\  nonlinearly  to  produce  secondary  sound 
waves.  The  secondary  wave  at  the  low  difference  frequency  wj  -  tog  passes  easily  through  the  relatively 
thin  tube  walls  and  radiates  to  the  far  field.  In  this  report  th't  far-fieid  >adiation  at  the  difference  frequency 
is  calculated  by  use  of  the  virtual-source-density  method.  Two  speciai  cases  are  discussed,  in  the  first  the 
high  frequencies  are  neighboring  resonance  frequencies  so  that  the  lengtl  of  SWAPS  is  nearly  1/2  wavelength 
at  the  difference  frequ  ncy  in  the  liquil  in  the  tube,  in  the  second  both  high  frequencies  lie  within  the 
bandwidth  of  a  single  .esonance  of  the  cavity,  in  this  second  case  SWAPS  can  be  as  short  as  1/4  wavelength 
at  the  high  frequency.  The  results  are  expressed  in  terms  of  a  quality  factor  that  describes  the  resonance 
properties  of  the  system  when  difference-frequency  radiation  is  oeing  generated.  This  quality  factor,  which 
must  be  determined  experimentally,  is  an  input  to  the  mathematical  me  del.  Comparison  is  made  to  a 
traveling-wave  acoustic  parametric  source  which  employs  the  same  piston  transducer  as  a  pump. 
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ABSTRACT 


A  standing-wave  acoustic  parametric  source  (SWAPS)  is  proposed  as  a  small, 
low-frequency  underwater  sound  source.  SWAPS  is  a  liquid-filled  cylindrical  tube 
which  is  driven  at  one  end  by  a  piston  transducer  and  terminated  on  the  other  end 
by  a  pressure-release  reflector  to  form  a  resonant  cavity.  The  piston  is  driven  simul¬ 
taneously  at  two  high  frequencies  and  o>2  which  are  at  or  near  resonance  for 
Pie  plane-wave  mode  in  the  cavity.  The  resulting  large-amplitude  primary  sound 
waves  mix  nonlinearly  to  produ  :e  secondary  sound  waves.  The  secondary  wave  at 
the  low  difference  frequency  oj1  —  u>2  passes  easily  through  the  relatively  thin  tube 
walls  and  radiates  to  the  far  field.  In  this  report  the  far-field  radiation  at  the  dif¬ 
ference  frequency  is  calculated  by  use  of  the  virtual -source-density  method.  Two 
special  cases  are  discussed:  in  the  first  the  high  frequencies  are  neighboring  reso- 
lance  frequencies  so  that  the  length  of  SWAPS  is  nearly  1/2  wavelength  at  the  dif- 
lt-ence  frequency  in  the  liquid  in  the  tube;  in  the  second  both  high  frequencies  lie 
within  the  bandwidth  of  a  single  resonance  of  the  cavity;  in  this  second  case  SWAPS 
can  be  as  short  as  1/4  wavelength  at  the  high  frequency.  The  results  are  expressed 
in  terms  of  a  quality  factor  that  describes  the  resonance  properties  of  the  system 
when  difference-frequency  radiation  is  being  generated.  This  quality  factor,  which 
must  be  determined  experimentally,  is  an  input  to  the  mathematical  model.  Com¬ 
parison  is  made  to  a  traveling-wave  acoustic  parametric  source  which  employs  the 
same  piston  transducer  as  a  pump. 
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MATHEMATICAL  MODEL  FOE  A  STANDING-WAVE  ACOUSTIC 
PARAMETRIC  SOURCE  (SWAPS) 


INTRODUCTION 

Within  the  past  decade  the  generation  of  narrow-beam,  low-frequency  acoustic  radia¬ 
tion  from  a  relatively  small  piston  source  has  been  accomplished  by  the  use  of  the  traveling- 
wave  acoustic  parametric  source.  Here  the  piston  source  is  driven  simultaneously  by  two 
primary  signals  of  high  frequency  and  co2.  During  subsequent  propagation  a  low- 
frequency  secondary  wave  at  the  difference  frequency  Wj  —  co2  is  generated  by  the  non¬ 
linear  interaction  of  the  primary  sound  waves  in  the  water  medium.  This  so-called  endfire 
array,  first,  proposed  by  Westervelt  [1]  and  confirmed  experimentally  by  Beilin  and  Beyer 

[2]  ,  has  been  the  subject  of  extensive  theoretical  and  experimental  investigation  by  Berktay 

[3]  and  others  [4-6].  The  desirability  of  such  a  device  is  enhanced  by  the  lack  of  side- 
lobes  in  the  radiation  pattern  and  by  the  relatively  large  oandwidlh  at  the  difference 
frequency.  Its  chief  disadvantage  is  the  low  efficiency  of  the  conversion  due  to  the  absorp¬ 
tion  and  diffraction  of  the  primary  waves  and  the  nonlinear  generation  and  subsequei  t 
absorption  of  both  the  netmonics  of  the  primary  waves  and  the  sum  frequency  tuj  +  ou2. 

This  limitation  is  unlike.y  to  be  eliminated  by  future  work. 

In  1965  Dunn,  Kuljio,  and  Welsby  [7]  briefly  examined  the  generation  jf  a  47-kHz 
sound  from  two  primary  waves  at  326  kFi  and  373  kHz  in  a  spherical-standing-wave 
system.  However,  they  were  primarily  concerned  with  subharmonic  generation  when 
cavitation  is  present  and  did  not  pursue  the  matter  further. 

In  this  report  we  describe  a  low-frequency  sound  source  called  a  standing-wave  acoustic 
parametric  source  (SWAPS)  which  increases  the  nonlinear  conversion  by  “folding”  the 
traveling  wave  to  produce  a  standing  wave  [8] .  This  is  accomplished  by  nse  of  a  closed 
cylindrical  tube  that  is  near  resonance  for  both  primary  sound  waves.  Because  the  parametric 
effect  is  inherently  nonlinear,  an  increased  generation  of  sound  at  the  difference  frequency 
results.  This  difference  frequency  passes  easily  through  the  relatively  thin  tube  walls  a’->d 
ends  and  radiates  to  the  far  field.  We  calculate  the  far-field  radiation  from  SWAPS  by  use 
of  the  virtual-source-densitj,  method  used  by  Westervelt  [1] .  Inc'uded  in  our  discussion 
is  a  comparison  to  a  traveling-wave  acoustic  parametric  source  which  employs  the  same 
piston  transducer  to  generate  the  primary  waves. 


DESCRIPTION  OF  SWAPS 

SWAPS  is  envisioned  to  be  a  liquid-filled  circular  cylindrical  tube  or  cavity  of  inner 
radius  a  and  lengtn  2  with  a  piston  transducer  mounted  on  one  end  and  a  pressure-release 
reflector  mounted  on  the  other  end  (Fig.  1).  The  piston  is  driven  simultaneously  at  the  two 
primary  frequencies  oo3  and  to2.  The  walls  are  made  sufficiently  rigid  to  allow  predominantly 
plane-wave  propagation  of  the  primaries.  Both  oj3  and  to2  are  chosen  to  be  near  a  resonance 
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Fig.  1  —  Standing-wave  acoustic  parametric  source 


frequency  of  the  plane-wave  mode  of  the  cavity.  Two  special  cases  art  considered:  (1) 
cu1  and  <jl>2  tiro  neighbor,  tg  resonance  frequencies;  (2)  uq  is  close  enough  to  co2  so  that 
they  both  he  well  within  the  bandwidth  of  a  single  resonance  frequency.  The  resulting 
primary  waves  present  in  this  standing-wave  system  now  interact  nonlinearly  to  produce 
sum  and  difference  components  as  well  as  various  harmonic  components.  The  difference 
frequency  corf  -  uq  —  u2  is  chosen  to  be  m  tch  smaller  than  either  primary  frequency. 

Thus  the  relatively  tlhn  walls  and  ends  of  the  tube  can  also  be  made  nearly  transparent 
acoustically  at  the  difference  frequency,  so  that  the  difference  frequency  radiates  unimpeded 
to  the  far  field. 

the  use  of  a  pres.su re-relense  reflector  should  inhibit  the  growth  of  both  the  sum 
frequency  and  harmonic",  of  the  primaries.  This  occurs  because  of  the  180°  phase  shift 
that  each  frequency  component  in  a  plane  wave  undergoes  upon  reflection  from  a  pressure- 
release  boundary.  Consider  a  sinusoidal  wave  leaving  the  piston.  As  it  propagates,  harmonics 
are  generated  with  a  fixed  phase  relationship  to  the  fundamental.  Ignoring  dispersion,  this 
phase  relationship  is  maintained  while  the  “most  stable”  waveform,  the  sawtooth,  is 
approached.  When  the  distorted  wave  is  reflected  from  the  pressure  rele~.se  end,  the  180° 
phase  shift  produces  a  “least  stable”  waveform  or  reverse  sawtooth  During  subsequent 
propagation  back  to  the  piston,  new  harmonic  generation  cancels  the  existing  harmonic 
convent,  and  the  waveform  ends  to  retun  to  a  sinusoid.  The  sum  frequency  cos  =  uq  + 
to2  will  be  close  in  frequency  to  the  second  harmonic  of  both  primaries  and  will  be  in¬ 
hibited  in  a  similar  manner.  In  addition,  any  dispersion  that  exists  in  the  tube  will  tend 
to  inhibit  the  growth  of  harmonics  and  the  sum  frequency,  since  it  will  also  tend  to  destroy 
the  stable  phase  relationship  between  the  various  components.  Thus  competing  nonlinear 
interactions  are  inhibited,  and  more  primary  energy  remains  available  for  generation  of  the 
difference  frequency. 

The  liquid  used  in  SWAPS  can  be  chosen  to  optimize  the  performance.  For  example 
an  ideal  liquid  m  ght  have  a  characteristic  acoustic  impedance  equal  *  j  that  of  the  surround- 
i  <g  water  medium  so  that  the  difference  frequency  is  transmitted  edactively  into  the  water, 
have  a  high  degree  of  nonlinearity  to  enhance  the  interaction  of  the  primaries,  have  small 
attenuatio  i  of  the  primaries  due  to  linear  loss  mechanisms  such  as  viscous  absorption  in 
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order  to  enhance  the  resonance  properties  o'  the  cavity,  have  a  high  cavitation  threshold 
(if  cavitation  proves  to  be  undesirable),  and  have  a  iow  sound  speed  to  maximize  the 
acoustic  size  of  the  interaction  volume. 


MATHEMATICAL  MODEL  OF  SWAPS 

The  basic  equations  of  acoustics  describing  the  conservation  of  mass,  momentum, 
and  energy  and  the  equation  of  state  of  the  medium  are  inherently  nonlinear.  Under  the 
assumptions  of  infinitesimal  wave  motion,  the  nonlinear  terms  can  be  neglected,  and  the 
linear  wave  equation  is  obtained.  Fortunately,  for  most  applications,  infinitesimal  theory 
is  quite  adequate.  However,  when  the  wave  motion  is  of  large  amplitude  so  that  the 
acoustic  particle  velocity  is  not  negligible  whan  compared  with  the  sound  speed,  the  non¬ 
linear  terms  must  be  retained.  is  the  existence  of  these  nonlinear  terms  that  leads  to 
mixing  of  the  primaiy  waves  to  produce  secondary  waves  which  consist  of  harmonics  of 
each  primary  frequency  component  and  intermodulation  components  such  as  sum  and 
difference  frequencies.  References  9  through  12  giv.  detailed  discussions  of  nonlinear 
acoustics. 

Because  of  the  intractable  nature  of  the  nonlinear  equations  one  usually  makes  a 
number  of  simplifying  assumptions  or  idealizations  to  obtain  a  solution.  For  example,  if 
one  assumes  that  the  secondary  waves  are  small  compared  with  the  primary  waves,  first- 
order  perturbation  theory  can  bo  used.  Such  an  approach  was  used  by  Westervelt  [1]  m 
describing  the  parametric  acoustic  array,  which  wc  call  a  traveling-wave  acoustic  parametric 
source  and  apply  the  acronym  TWAPS. 

Westervelt  used  as  his  starting  point  the  exact  equation  for  arbitrary  fluid  motion 
derived  by  Lighthill  [13] .  Using  perturbation  theory  and  neglecting  viscosity,  he  derived 
the  inhomogeneous  wave  equation 


Q2p  =  -/>o  (i) 

where  P  is  the  acoustic  pressure  of  the  secondary  wave  which  includes  the  harmonic  and 
intermodulation  components  and  p0  *s  ^ie  equilibrium  mass  density.  The  virtual-simple- 
source-strength  density  q  which  results  from  the  primary-wave  pressure  Pi  is  given  by 


<1  = 


p2  3t 
fJQc0 


(Pf), 


(2) 


where  is  the  sound  speed  for  infinitesimal  waves.  The  parameter  T  is  defined  by 

P  =  i  +  po_ 

2c2  \dp2  )f>=P0 


f3) 


or 


r  = 


(4) 
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where  the  nonlinearity  parameter  B/A  is  a  measure  of  the  nonlinear  response  of  the  liquid 
medium.  For  water  at  20°,  B/A  is  about  5.  Other  liquids  have  values  of  B/A  that  lie 
b  :  tween  about  4  and  11  [11] . 

Westervelt  assumed  a  primary  wave  consisting  of  two  collimated  collinear  plane  waves 
of  frequencies  coj  and  co2  emanating  from  a  common  source.  He  reintroduced  the  effects 
of  viscosity  in  an  ad  hoc  way  by  assuming  that  the  primaries  are  attenuated  by  the  usual 
absorption  coefficient  of  linear  acoustics.  Substituting  this  primary  wave  into  Eq.  (2)  led 
to  a  virtual-source  density  with  the  frequency  components  2co1}  2u>2,  o)s,  and  ood.  He 
then  obtained  the  far-field  radiation  at  the  difference  frequency  by  using  the  cod  virtual- 
source-density  term  to  obtain  the  inhomogeneous  term  f  in  Eq.  (1)  and  by  integrating 
this  expression  times  the  free  space  Green’s  function  over  the  cylindrical  volume  containing 
the  primary  waves.  He  assumed  that  the  difference  frequency  is  low  enough  so  that  its 
attenuation  is  negligible.  This  model  assumes  that  the  transfer  of  energy  from  the  primaries 
into  secondary  waves  such  as  harmonics  and  sum  and  difference  components  is  small 
enough  so  that  nonlinear  attenuation  of  the  primaries  is  negligible.  Recently  attempts 
have  been  made  to  modify  the  primary-wave  distribution  to  approximate  the  nonlinear 
attenuation  due  to  harmonic  generation  [6] .  At  large  primary  amplitudes  (or  at  very  low 
primary  frequencies),  this  loss  mechanism  will  dominate  linear  absorption. 

We  now  calculate  the  difference  frequency  radiated  from  SWAPS  into  the  far  field 
by  the  perturuation  approach  of  Westervelt.  First  we  obtain  the  primary -wave  distribution. 
The  piston  end  of  SWAPS  is  assumed  to  be  vibrating  with  a  normal  velocity  V  given  by 

V  =  V1  sin  cujt  +  H2  sin  co2t.  (5) 

The  axis  of  the  tube  is  taken  to  be  the  z  axis>  with  the  piston  at  z  =  0.  The  reflector  end 
of  the  tube  is  pressure  released,  so  that 

P{  =  0  at  2  =  Si.  (6) 

The  walls  of  the  tube  are  assumed  rigid  at  Wj  and  to2,  and  plane  wave  motion  is  assumed. 
The  linear  attenuation  coefficients  for  the  primaries  in  the  tube  ctl  and  a2  include  contri¬ 
butions  from  such  linear  loss  mechanisms  as  boundary-layer  and  mainstream-viscous 
attenuation.  The  solution  to  the  one-dimensional  linear  wave  equation  subject  to  the 
boundary  conditions  given  in  Eqs.  (G)  and  (6)  is  given  by 


=  p0  V1  Re 


(co1  sin  [(£  —  z)  (fe1  —  ) 

|  (fti  -  telcos  [2(fej  -toj)]  j 


+  a  similar  term  for  P2 , 


(7) 


where  Re  indicates  the  real  part  of  the  expression  in  braces  and  k  is  the  wavenumber. 


We  now  assume  that  the  linear  attenuation  is  small,  so  that 


ajK  «  1, 


(8) 
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and 


a2£«l.  (8) 

In  addition  we  assume  that  the  length  £  is  such  that  both  Wj  and  oj2  are  equal  to  or 
nearly  equal  to  a  resonance  frequency  of  the  standing-wave  system,  that  is, 

OJj  —  0Jn  «  con 


and 


(9) 


“>n  ~  “>2  « 

where  it  is  assumed  that  and  co2  are  located  on  opposite  sides  of  con.  When  the  small 
dispersion  and  attenuation  that  will  exist  in  the  tube  are  neglected,  the  resonance  fre¬ 
quencies  con  are  given  by 


_  (2n  +  1)tt Cq 
“n  2£ 


(10) 


where  the  integer  n  indicates  the  order  of  the  resonance  and  c0  is  the  free-field  sound 
speed  for  the  liquid  in  SWAPS. 

With  these  assumptions  the  linear  solution  can  be  written  in  the  simplified  form 
Pj  -  ,o0Cq  [ V1Al  sin  {(Ojf  —  0t)  cos  kxz  +  V2A2  sin  (c o2t  —  02)  cos  k2z]  >  (H) 
where  the  amplitude  and  phase  factors  are  given  by 


CVjS  (1  +  [(Wj  -  )/(£*! c0)] 2)  i/2 


and 


(12) 


01  =  tan  -  w„)/(a1c0)], 

and  Similarly  for  A2  and  02-  At  resonance,  Wj  =  con  and  A 2  «ud  A2  have  their  maximum 
values  of  l/a^  and  l/a22  respectively.  From  Eqs.  (12)  we  see  that  the  conditions  that 
cjj  and  co2  lie  within  the  bandwidth  of  the  resonance  at  uin  can  be  expressed  as 


<*1  wn 

alc0 


<  1 


un  -  ^2 
a2c0 


<  1. 


and 


(13) 


6 


ROGERS  AND  VAN  BUREW 


If  the  amplitudes  of  the  primary  waves  were  small  and  if  the  assumed  boundary 
conditions  at  the  walls  and  ends  were  met,  Eq.  (11)  would  be  a  good  representation  of 
the  primary  waves.  However  these  primar  r  waves  would  not  produce  much  difference- 
frequency  radiation.  Therefore,  if  signific.  nt  difference-frequency  radiation  is  desired, 
then  the  amplitudes  of  the  primaries  must  be  large  and  the  infinitesimal  theory  used  to 
obtain  Eq.  (11)  is  no  longer  valid.  Even  i~  the  pressure-release  reflector  works  perfectly, 
the  harmonics  of  the  primaries  as  well  as  other  high-frequency  intermodulation  components 
cannot  now  be  treated  as  small  secondary  waves  which  have  a  negligible  effect  on  the 
primary  waves  and  can  be  calculated  from  first-order  perturbation  theory.  Thus  we 
modify  the  primary-wave  distribution  to  include  these  high-frequency  components  and  still 
consider  the  difference  frequency  cjd  to  be  the  secondary  wave.  Since  wo  are  interested 
in  only  u>d,  we  need  to  know  only  the  distribution  of  the  fundamental  components  of 
the  primaries.  We  assume  that  this  distribution  is  given  by 

Pj  =  PoCof^Qi  sin  (coxZ  —  <px)  cos  kxz  +  V2Q 2  sin  (co2t  —  02)  cos  k2z],  (14) 

where  the  quality  factors  Qt  and  Q2  contain  the  dependence  given  by  Eqs.  (12)  and  an 
additional  unspecified  factor  representing  the  loss  of  energy  from  the  fundamental  com¬ 
ponents  into  the  generation  of  high-frequency  intermodulation  components.  We  also  include 
in  Qx  and  Q2  t*ie  l°ss  energy  from  co1  and  to2  due  to  generation  of  the  difference 
frequency  itself.  This  apparent  departure  from  second-order  perturbation  theory  is 
justified  by  the  fact  that  the  difference  frequency  does  not  build  up  in  the  standing-wave 
system  but  instead  radiates  to  the  far  field.  Consequently  the  amount  of  difference  fre¬ 
quency  present  in  the  region  occupied  by  the  primaries  can  be  small  compared  with  the 
primaries  even  if  a  significant  amount  of  difference  frequency  is  radiated.  Therefore  Qj 
and  Q2  determine  the  amplitudes  of  the  primaries  during  difference-frequency  generation 
and  are  each  a  function  of  both  piston  velocity  components  Vl  and  V2.  To  simplify  the 
analysis,  we  assume  that  Vj  is  equal  to  V2  and  that  —  o>„  =  —  u>2  -  cod/2  when 

both  primaries  lie  within  a  single  resonance.  Then  since  co2  =  to,  so  that  a2  — 

we  have  rij  ^  A 2,  Q1  ^  —  Q,  and  0X  =  ~02  —  0o-  The  Primary  distribution  becomes 

Pi  ~  PoeoVGIsin  “  $0 )  cos  k\z  +  sin  (co2<  +  <p0)  cos  k2z] .  (15) 

Thus  we  require  that  the  behavior  of  Q  as  a  function  of  V  be  determined  experimentally 
as  an  input  to  this  mathematical  model. 

Substituting  the  expression  for  Pi  into  Eqs.  (1)  and  (2)  and  retaining  only  the  terms 
that  v  ontribute  to  the  difference  frequency,  we  have 


f-„  rp0v*<}2*3 

'  p°  dt~  4 


(cos  (c odt  —  kdz  —  20o )  +  cos  (c odt  +  kdz  —  20o ) 


+  cos  (c odt  —  ksz  —  200 )  +  cos  (c odt  +  ksz  —  20o)] . 


(16) 


It  can  be  shown  that  the  last  two  terms  do  not  contribute  significantly  to  the  far  field ^ 
when  cj,  »  oj({  and  ks  »  kd.  The  difference-frequency  pressure  at  the  field  point  R 
is  given  by 
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f-ikd  \R  -  ?\ 

|K  -r I 


dV, 


(17) 


where  t]ie  integration  is  over  the  circular  cylindrical  volume  of  radius  a  and  length  £  and 
where  r  is  the  source  point  (Fig.  2).  We  assume  for  thej>resent  that  SWAPS  is  filled  with 
water  (or  whatever  the  surrounding  liquid  is).  We  take  R  to  be  in  the  far  field  and  use 
the  asymptotic  form  for  the  free-space  Green’s  function.  The  tree-space  Green’s  function 
is  appropriate  for  the  difference  frequency,  since  we  assumed  that  the  walls  and  ends  of 
SWAPS  are  acoustically  transparent  at  cod.  Neglecting  the  last  two  terms  in  Eq.  (16), 
converting  to  complex  notation,  and  choosing  cylindrical  coordinates  (o,  j3,  z)  for  the 
source  point  and  spherical  coordinates  (R,  0,  <j>)  for  the  field  point,  we  obtain  the  expression 


ltdR  2 0O) 

Pd(R,0)  =SF(0)e- - ^ - ,  (18) 

with 

SF(0)  =-•  f  (Q  eikAo^0cos(^3)^os0}{~ikdZ  +  eikdz)  odadpdZt 

16?r  Jo  Jo  Jo 

(19) 


where  S  is  the  source  level  and  F(0)  is  .he  far-field  pressure  distribution.  The  integration 
over  /?  of  the  factor  depending  on  0  and  a  produces  27rJo(/cda  sin  0 ),  where  is  the 
zero-order  Bessel  function,  and  subsequent  integration  over  u  produces  the  aperture  factor 
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H  =  kfito  Ji  ( V sin  «>•  <2°) 

We  usually  choose  kda  sin  0  «  1,  so  that  J1  can  be  replaced  by  its  lowest  order  term 
and  H  has  its  maximum  value 


H  ~  7 ra2. 


(21) 


Integration  over  z  of  the  remaining  factors  produces  two  oppositely  directed  end-fired- 
array  terms  that  are  kd2  out  of  phase  with  each  other: 


E<0)  =  eik ‘i*  f  Sm  C°s2  +  fsin  IM  sin2(0/2)]  ] 

{  /erf2  cos2(0/2)  j  \  /erf2  sin2(d/2)  j 


(22) 


The  second  term  in  braces  can  be  obtained  from  the  first  term  in  braces  by  replacing  0 
with  180°  —  0.  Collecting  terms,  we  obtain  the  source  level 

S  =  rp0V-Q2^tt22|£(0)l/16  (23) 

anC  the  far-field  pressure  distribution 


J-,  (k,i a  sin  0) 

|F(0)|  =  - t-4  |£(0)|- 

'  {l/2)Kda  sm  0  v  n 


sin2[/erf2  cos2(0/2)]  sin2[/?d2 


am 


+  2  COS  k([2 


cos'1  (d/2)  sin4  (d/2) 

1/2 


sin  [kd2  sin2(d/2)j  sin  [7jj2cos2(0/2)] 
sin2(0/2)  cos2(0/2) 


(24) 


where 


iE(0)  I  = 


sin2/erf£ 

ft2  2  2 
d 


sin  2/e,i2 

+ - 

kdl 


(25) 


and  where  |F(0)|  is  normalized  to  unity  in  the  forward  direction  (0  =  O' 

m  the  first  special  case,  where  and  cu2  are  two  neighboring  resonance  frequencies, 
it  can  easily  be  shown  that  kd2  =  it,  so  that  2  =  \dl 2.  Then  !F(0)|  is  equal  to  unity,  and 
the  source  level  and  fa'-field  pressure  distribution  for  this  half-wavelength  SWAPS  (HWS) 
are  given  by 


Shivs  =  7rPp0  V2Q2/erfa2/16 


(26) 


and 


|F(0)|  = 


<.Ji(kda  sin  0) 

cos  0  cos  (7T/2  cos  0 ) 

(1/2  )kda  sin  0 

(7r/4)  sin20 

(27) 
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We  assume  that  the  first  factor  is  nearly  equal  to  unity  and  plot  the  second  factor  in 
Fig.  3.  The  far-field  pressure  distribution  becomes  zero  in  the  broadside  direction,  similar 
to  the  far-field  pressure  distribution  of  a  linear  dipole.  The  3-dB-down  points  occur  at 
0  =  53°  and  0  =  127°.  Complete  nulls  occur  in  this  pattern  because  the  two  oppositely 
directed  waves  in  the  virtual-source  density  produce  complete  cancellation  over  a  half 
wavelength.  If  SWAPS  contained  a  liquid  with  a  sound  speed  less  than  that  of  water, 
then  the  length  of  the  tube  would  be  less  than  a  half  wavelength  in  water  and  the  cancella¬ 
tion  would  be  incomplete.  This  liquid  should  also  have  a  density  that  is  correspondingly 
larger  than  that  of  water,  so  that  the  ebara  ;teristic  impedances  would  be  nearly  equal  and 
the  effect  of  the  boundary  between  the  two  liquids  would  be  minimized.  This  impedance 
match  coupled  with  the  small  acoustical  size  of  SWAPS  at  the  difference  frequency  should 
allow  the  far  field  to  be  calculated  as  if  SWAPS  were  filled  with  water  insofar  as  the  phase 
terms  are  concerned.  Of  course  the  quantities  v .  c0,  and  F  appearing  in  the  amplitude 
of  f  should  be  evaluated  for  the  actual  liquid  in  SWAPS.  For  example,  if  the  liquid  in 
SWAPS  were  similar  to  carbon  tetrachloride  with  a  sound  speed  of  105  centimeter  per 
second  and  a  density  of  1.5  grains  per  cubic  centimeter,  the  far-field  pressure  distribution 
calculated  using  Eq.  (24)  would  be  nearly  omnidirectional  with  only  1/2  dB  difference 
between  the  maxima  at  0°  and  180°  and  the  minima  at  90°  md  270°. 

The  total  power  W//  \vs  radiated  at  the  difference  frequency  by  the  half-wavelength 
SWAPS  is  obtained  by  numerical  integration  of  |F(0)|2  with  Gaussian  quadrature.  The 
result  is 


WHWS  ~ 


1-64  rrSfnvs 
2Poco 


(28) 


Fig.  3  —  Far-field  pressure  distribution  for 
the  half-wavelength  SWAPS 


10 


ROGERS  AND  VAN  BUREN 


Of  perhaps  more  interest  is  the  second  special  case,  where  both  cox  and  co2  lie  witl  in 
the  bandwidth  of  a  single  resonance  of  the  standing- wave  system.  In  this  case  the  length 
of  SWAPS  is  nearly  equal  to 


2  = 


X, 


(251) 


where  m  is  an  integer  and  X  is  the  wavelength  in  SWAPS  at  either  primary  frequency. 
Consequently  £  can  be  considerably  shorter  than  a  half  wavelength  at  the  difference 
frequency  to  take  advantage  of  the  fact  that  S  depends  on  Q22.  Thus  if  Q2  increases 
faster  than  1/2  as  2  is  decreased  to  values  "srresponding  to  smaller  integer  values  for  in, 
there  will  bt  an  overall  increase  in  S.  For  small-amplitude  waves  A  is  seen  from  Eq.  (12) 
that  Q 2,2  inc 'eases  nearly  as  1/2.  For  large-amplitude  waves,  Q°2  should  also  increase  as 
the  length  is  decreased  from  that  of  a  half-wavelength  SWAPS.  At  some  optimum  length, 
however,  increased  nonlinear  losses  possible  cavitation,  and  heating  effects  will  cause  Q2 2 
to  begin  to  decrease. 

For  th!s  short  SWAPS  (SS),  the  far-held  pressure  distribution  predicted  by  Eq.  (24)  is 
omnidirectional,  and  the  source  level  and  radiated  power  are  given  by 

Sss  =  Vp0V2Q2k2da2m,  (00) 


and 


W'SS 


_  ^Sss 

2 PqCq' 


(31) 


To  get  some  feeling  for  these  expressions,  we  compare  the  results  for  the  short  SWAPS 
to  those  obtained  by  Westervelt  for  TWAPS,  although  we  recognize  that  SWAPS  and  TWAPS 
are  designed  for  different  applications.  For  simplicity  we  assume  that  TV/APS  is  operating 
with  the  same  piston  velocity  as  SWAPS.  The  source  level  and  radiated  power  calculated 
by  Westervelt  are  designated  by  .S-/>  and  WT,  and  w<  use  o'  to  represent  the  free-field 
attenuation  of  the  primaries  in  seawater.  This  attenuation  coefficient  should  include  con¬ 
tributions  from  viscous  attenuation,  the  generation  of  harmonics  and  other  intermoduiatum 
components,  and  diffraction.  For  the  case  where  the  short  SWAPS  is  filled  with  water, 
we  have 

"  32«'2  (32) 

Orp 

and 


~  =  2Q4a'2Arf2/7r.  (33) 

WT 

Thus  the  achievement  of  a  source  level  comparable  to  that  of  TWAPS  requires  that  Q 2  be 
comparable  to  l/a'2.  We  are  nut  predicting  values  for  Q  at  the  present  time  and  emphasize 
again  that  it  is  a  parameter  that  must  be  experimentally  determined  as  an  input  to  this 
mathematical  model.  The  ratio  of  total  radiated  powers  is  much  more  favorable  than  the 
ratio  of  source  levels,  since  kd  is  expected  to  be  much  larger  than  a\  so  that 
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^SS 

*1 


» 


(33) 


SUMMARY 

We  summarize  some  of  the  more  important  features  of  the  standing-wove  acoustic 
parametric  source  (SWAPS): 

•  SWAPS  is  designed  to  be  a  small,  low-frequency  underwater  sound  source. 

•  Two  special  modes  of  operation  are  available: 

-The  two  primary  frequencies  are  neighboring  resonance  frequencies  for  the 
standing-wave  system  in  SWAPS.  In  this  case  the  lengtn  of  SWAPS  is  a  half  wavelength 
at  the  difference  frequency  in  the  liquid  in  SWAPS. 

—The  two  primary  frequencies  are  both  within  the  bandwidth  of  a  single 
resonance  frequency  for  the  standing-wave  system  in  SWAPS.  In  this  mode  the  length 
of  SWAPS  is  nearly  equal  to  an  odd  number  of  quarter  wavelengths  at  either  primary 
frequency  in  the  liquid  in  SWAPS  and  thus  can  be  considerably  shorter  than  the  half¬ 
wavelength  SWAPS. 

®  SWAPS  is  expected  to  be  a  more  effective  radiator  of  difference  frequency  than 
a  traveling  wave  acoustic  parametric  source  (TWAPS)  because  (a)  high-energy  density  and, 
consequently,  increased  nonlinear  conversion  is  achieved  by  simultaneous  resonance  or 
near  resonance  of  the  two  primary  signals  and  (b)  competing  nonlinear  interactions  such 
as  the  generation  of  harmonics  of  the  primaries  and  the  generation  of  the  sum  frequency 
are  inhibited  by  the  pressure-release  reflector. 

»  Since  SWAPS  is  short  acoustically  at  the  difference  frequency,  it  will  not  have  the 
high  directivity  associated  with  TWAPS.  Instead  SWAPS  is  intended  for  applications  that 
require  an  omnidirectional  or  perhaps  a  dipolar  far-field  pressure  distribution.  If  directivity 
is  desired,  an  array  of  SWAPS  sources  can  be  used.  The  frequently  encountered  acoustic 
interaction  problem  inherent  with  an  array  of  small  conventional  sources  should  not  occur 
here  because  the  piston  which  produces  the  high-frequency  primary  waves  in  a  SWAPS 
source  will  not  be  affected  appreciably  by  the  low-frequency  difference-frequency  radiation 
from  neighboring  SWAPS  sources  The  high-frequency  primary  waves  are  contained  within 
the  SWAPS  sources  and  therefore  cannot  affect  the  piston  in  a  neighboring  SWAPS  source. 

•  The  liquid  used  in  SWAPS  can  be  chosen  to  optimize  the  performance.  Some 
important  parameters  to  consider  are  sound  speed,  density,  nonlinearity,  attenuation 
coefficient,  and  cavitation  strength. 


ACKNOWLEDGMENT 

The  authors  gratefully  acknowledge  the  continued  support  of  NavShips:  PMS302-4. 


12 


ROGERS  AND  VAN  BUREJJ 


REFERENCES 

1.  P.J.  Westervelt,  “Parametric  Acoustic  Array,”  J.  Acoust.  Soc.  Amer.  35,  535-537  (1963). 

2.  J.L.S.  Beilin,  and  R.T.  Beyer,  “Experimental  Investigation  of  an  End-Fire  Array,”  J. 
Acoust.  Soc,  Amer.  34,  1051-1054  (1962). 

3.  H.O.  Berktay,  “Possible  Exploitation  of  Non-Linear  Acoustics  in  Underwater  Trans¬ 
mitting  Applications,”  J.  Sound  Vib.  2,  435-461  (1965). 

4.  R.H.  Mellen,  W.L.  Konrad,  and  D.G.  Crowning,  “Approximate  Scaling  Laws  for 
Parametric  Sonar  Transmitter  Design,”  Proceedings  of  the  British  Acoustical  Society 
Meeting  on  Nonlinear  Acoustics,  Birmingham,  England,  Apt.  1971. 

:  5.  T.G.  Muir,  “An  Analysis  of  the  Parametric  Acoustic  Array  for  Spherical  Wave  Fields,” 
Applied  Research  Lab.,  Univ.  Texas,  Tech.  Report  ARL-TR-71-1,  May  1971. 

6.  R.H.  Mellen,  and  M.B.  Moffett,  “A  Model  for  Parametric  Sonar  Radiator  Design,” 

U.S.  Navy  Journal  of  Underwater  Acoustics  22,  105-116  (1972). 

7.  D.J.  Dunn,  M.  Kuljis,  and  V.G.  Welsby,  “Non-linear  Effects  in  a  Focused  Underwater 
Standing  Wave  Acoustic  System,”  J.  Sound  Vib.  2,  471-476  (1965). 

8.  A.L.  Van  Buren,  and  P.H.  Rogers,  “Mathematical  Model  for  a  Standing  Wave  Acoustic 
Parametric  Source,”  J.  Acoust.  Soc.  Amer.  53,  382  (A)  (1973). 

9.  R.T.  Beyer,  “Nonlinear  Acoustics,”  chapter  10  in  Physical  Acoustics,  Vol.  II,  Part  B, 
edited  by  W.P.  Mason,  Academic  Press,  New  York,  1965. 

10.  D.T.  Blackstock,  “Nonlinear  Acoustics  (Theoretical),”  chapter  3n  in  American  Institute 
of  Physics  Handbook,  3rd  edition,  edited  by  D.E.  Gray,  McGraw-Hill,  New  York,  1972. 

11.  R.T.  Beyer,  “Nonlinear  Acoustics  (Experimental),”  chapter  3o  in  American  Institute 

of  Physics  Handbook,  3rd  edition,  edited  by  D.E.  Gray,  McGraw-Hill,  New  York,  1972. 

12.  L.D.  Rozenberg,  High-Intensity  Ultrasonic  Fields,  Plenum  Press,  New  York,  1971. 

13.  M.J.  Ligh thill,  Math.  Revs.  19,  915  (1958). 


